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We consider two different kinds of fluctuations in an ion trap potentiai: external fluctuating 
electrical fields, which cause statistical movement ("wobbling") of the ion relative to the center of 
the trap, and fluctuations of the spring constant, which are due to fluctuations of the ac-component 
of the potential applied in the Paul trap for ions. We write down master equations for both cases 
and, averaging out the noise, obtain expressions for the heating of the ion. We compare our results to 
previous results for far-off resonance optical traps and heating in ion traps. The effect of fluctuating 
external electrical fields for a quantum gate operation (controlled-NOT) is determined and the 
QQ ■ fidelity for that operation derived. 
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I. INTRODUCTION 



The ability to engineer and control pure quantum states of trapped ions is driving a number of new technologies 
including time and frequency measurements, new measurement schemes and quantum logic manipulations for quantum 
computation Q-^]. The key requirement for engineering quantum states is the necessity to remove, or at least control, 
all sources of noise and uncertainty. Laser cooling in particular enables uncertainty about the initial vibrational state 
to be removed by cooling the ions to the collective vibrational ground state. From that point pure states may be 
prepared using highly stabilized laser pulses. Despite these achievements, however, noise cannot be entirely eliminated. 
Residual laser intensity and phase fluctuations in the pulses that are used to shape the quantum states must be taken 
into account as well as noise in the trapping parameters. In this paper we consider various sources of noise in 
the trap itself and determine their effect on cold trapped ions and the ability to perform reversible logical operations. 
. Noise is of course the origin of decoherence, the process which limits the ability to maintain pure quantum states. 
However, we need to be careful in making the connection between noise and decoherence. From a fundamental 
perspective the dynamics of the ions is always unitary and reversible, even in the presence of noise, but by definition 
the dynamics of a noisy quantity is uncontrollable and often unknown. Thus the precise unitary dynamics varies 
from one run of the experiment to the next and the exact motion of the state in Hilbert space may not be known or 
even precisely predictable. In the course of the experiment we do not have precise control over, or knowledge of, the 
unitary transformations in state engineering, and thus we cannot be sure we have reached the desired state in Hilbert 
O^' space. In the case of quantum computation such a result is manifest as an error. Faced with describing such a system 
we can simply average over the noise, which in practical terms means we combine the data from many experiments 
all performed with different realizations of the noisy control parameters. Alternatively we can give the sample space 
of error states in each run, together with their probability of occurrence. In this paper we give both descriptions 
particularly for the case of fluctuations of the center (equilibrium point) of the trap potential. 

In an ion trap an inability to precisely control the motion leads to an unwanted excitation of the vibrational state 
of the ion, that is 'heating'. The main source of this heating appears to be due to the ambient fluctuating electrical 
fields in the trap. There are now a number of experiments that have measured this heating Recently James 

H has shown that a simple theory of this source of heating can be given in terms of a harmonic oscillator subject to a 
fluctuating classical driving field. Our model for this source of heating is similar, although our theoretical description 
is a little different. 

The paper is organized as follows: In section II we first derive the heating rate of the ion due to fluctuating electrical 
fields. Section III is devoted to the effects of fluctuations in the spring constant of the trap potential on the heating. 
In the fourth section we look at the effects of fluctuating electrical fields on gate operations. As a specific example we 
investigate the effects on the so far experimentally realized controlled-NOT gate (NIST gate) pl| . We conclude with 
a discussion of our results. 



II. HEATING DUE TO FLUCTUATING ELECTRICAL FIELDS 

We want to model the effects of fluctuating electrical fields with the same formalism used in [Q and compare our 
results to those obtained by Savard et al. (h]] . There they derive the heating rate for a far-off resonance optical trap 
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due to fluctuations in the location of the trap center which are caused by laser-beam-pointing noise. We apply the 
formalism here to an ion in a Paul trap. In those traps additional electrical fields cause a replacement of the center 
of the trap, to which the ion adjusts automatically by finding the minimum of the potential. We assume that those 
electrical fields which cause this replacement have got an additional white noise component, to which the ion cannot 
adjust. 

The Hamiltonian in this case is 

H = £ + \ muj 2x2 + mx ' (1) 

where x and p are the position and momentum operator of the ion in the trap, m is the mass of the ion, and u) is the 
trap frequency. The term describes the fluctuating electrical force due to fluctuating fields which we assume to 
be due to a white noise process, i.e. 

C(t)dt = VldW(t) , (2) 

with dW(t) a Wiener process and the parameter 7 scales the noise. Since we are dealing with a white noise process, 
we define a stochastic Schrodinger equation in the Ito formalism [n2| 



dp = -^[H ,p\dt-^[ X ,p}dW(t)-^[x,[x,p\] , (3) 



where 



H = Z- + \rruJx 2 . (4) 

For a single run with a known (or monitored) behaviour of the noise in time, the above equation gives the evolution 
of the system density operator conditioned on the entire history of the noise process. Since we are not interested in 
the effects of the fluctuations in some short time limit or just for one run of the experiment, we follow H and average 
over the noise to get the master equation for the average density operator 



dp »,„ 7 



- T [H ,p\ - ^2 [at, [x,p\] . (5) 



This equation has a 'high frequency' limit which may be relevant in experiments for which the time scale is much 
longer than the period of the trap. For example if the trapped ions form a quantum logic gate the time over which 
the gate operation is imposed may be much greater than the trap period (which is typically of the order of 10 _6 s). 
In that case we can transform to a frame rotating at the trap frequency and time average the rapidly rotating terms 
to give a master equation in the form 

dp 7 / + t 1 / t t t t\\ 

apa' + a' pa — —(a ' ap + aa 'p + pa 'a + paa ') I . (6) 



dt 2hmui V 2 
To calculate the mean energy 



{H ){t) = ^ + \mu 2 {x 2 ) , (7) 
2m 2 

we look at the time derivative of (x 2 ) and (p 2 ) and get the two equations 

(x 2 ) = — — (xp + px) (8) 

TO 

(p 2 ) = muu 2 {xp + px) + 7 , (9) 

and thus 

d(H ) _ 1 



7 ■ (10) 
dt 2m' y ' 



So the mean energy is 
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(ff )(t) = -! 7 t + <if ) (t=0) . (ii) 

Ira 

(The same result obtains with the time averaged master equation as well). This result is equivalent to that derived 
by Savard et al. fiof and by James Q in the limit of white noise and with appropriate changes of notation. 

To make the comparison with the results of Wineland et al. Q we write Eq. (|ll]) in terms of the mean vibrational 
quantum number n as 

where the time constant for decoherence is then given by 

2hmuj , 
t* = . (13) 

7 

The fluctuating linear potentials are caused by fluctuating electric fields, E(t), on the trap electrodes thus we expect 
that the fluctuating term in Eq. (|l) is given by £(t) = qE(t). If the fluctuations in the electric field are treated as 
white noise the spectral density of these fluctuations (near the trap frequency) is independent of frequency. If we take 
E(t)dt = EodW(t) then we have the equivalence 7 = q 2 E 2 with the spectral density of fluctuations in the field given 
by 

S e {lo) = 4 / E(t + T)E{t)dr (14) 
Jo 

= 2£ 2 • (15) 



Thus the decoherence time becomes 



which is the form quoted in reference g] . 
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III. NOISE IN THE SPRING CONSTANT 



In an ion trap a suitable combination of ac and dc electric fields are used to create an approximate harmonic 
potential in three spatial dimensions for one or more ions |l7| ]. In the case of more than one ion, the coulomb force 
couples the motion of the ions and the collective normal mode coordinates undergo harmonic motion. Laser cooling 
enables the ion, or the collective mode of many ions, to be prepared in or near the ground state of the system. Noise 
in the spring constant in due to fluctuations in the ac-component and dc-components of the applied potential. In the 
linear trap of reference [^) with the long axis oriented along the z-axis, and with micromotion ignored, the periodic 
trap potential at the trap centre is given by 

$ « -rm? z + ^mcj 2 (x 2 + y 2 ) , (17) 

where the harmonic force in the z-direction is formed by a static potential and in the x — y direction it is formed by an 
ac potential at frequency £It and amplitude Vq. The resonance frequency is then given by uj r = qV /{2 1 / 2 n T mR 2 ), 
where m and q are the ion mas and charge and R is the distance from the z-axis to the surface of the linear electrodes. 
Clearly fluctuations in either the dc or ac component will lead to fluctuations in the spring constants of the trap, 
although in practice ac fluctuations are more significant. 

In far-off resonance optical traps, the restoring force is provoced by the induced optical dipole force of an applied 
laser. Typically the atom is confined at the node of a standing wave of a laser tuned above the atomic resonance 
(blue detuning). The atom them sees a mechanical potential proportional to the intensity of the laser. If the intensity 
is quadratic near the node, a linear restoring force will be produced. In this case fluctuations in the applied laser 
intensity lead to fluctuations in the trap frequency | |L0| . 

The Hamiltonian for a particle moving in a harmonic potential with fluctuating spring constant is 

H=^ + \ m uj 2 {l + e{t))x 2 , (18) 
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where again we assume white noise, 

e(t)dt = VTdW(t) , (19) 
and the noise this time is scaled with the parameter T. To make things easier, we introduce dimensionless coordinates 

-1/2 



x = X 



\moj ) 
p — > (2hmuj)^ 1 ^ 2 p = P , 



(20) 
(21) 



so that the commutation relation in these new coordinates now reads 



Thus our Hamiltonian takes on the form 



[x,P} = - 2 . 



H = fiuj (P 2 + X 2 ) + e{t)huoX 2 
= H + e{t)TiujX 2 . 



(22) 



(23) 



Again we are not interested in a specific noise history, but rather in the overall effect of the fluctuations. So we 
average out the noise and the corresponding master equation reads 



(24) 



We want to determine the change of the mean energy with time. The mean energy is given by the expression 

(H a )(t) = hcj((P 2 )(t) + (X 2 )(t)) . (25) 

We wish to get the expressions for (P 2 }(t) and (X 2 )(t). To do this, we can first derive a system of first order differential 
equations for (P 2 )(t), (X 2 )(t) and (XP + PX)(t)/2: 



d 
~dt 



(X 2 ) 
(P 2 ) 
\{XP + PX) 



= A 



(X 2 ) 

(P 2 ) 

\{XP + PX) 



(26) 



where 



2uj 
Too 2 -2lu 
-uj uj 



(27) 



So we have to solve this system of differential equations to get the solution for (X 2 )(t) and (P 2 )(t). The exact solution 
is 



(H )(t)=hu< 



,2{D 2 - l)w (2 — D 2 + 2D 4 )(1 + D 2 + L> 4 ) 
exp( -j= 1)- 



+ exp(- 



V3D 
(D 2 - l)u 
V3D 



«)( 



9D 2 (1-D 2 + D i ) 



9D 2 (1 - D 2 + L> 4 ) 



D 

2(1 -D 2 )(l + 2D 2 -Z? 4 + L> 6 + D 8 ) , (1 + D 2 )lu 
H t^tz^tz sin(- 



exp( 
+ exp(- 
+ 



3V3L> 2 (1 + L» 6 ) 
2(D 2 - l)w 1 + D 2 + D 4 



(X 2 )t=o 



V3D ' t h{l-D 2 + D i ) 
(D 2 - l)w / 2{D 2 - I) 2 



V3Z? " 
2D 2 {D 2 - 1) 



3(l-D 2 + D 4 



cos(- 



(1 + D 2 )u 



D 



t) 



V3(l + D 2 )(l + D 2 + D 4 ) 



. . (i + d 2 )uj : 

sin(- D ' t) 



(P 2 )t= 



4 



, (D 2 - l)w 
+ exp(— - 



2{D 6 



3V3-D(1 — D 2 + L> 4 ) 
2(D 6 - 1) 
V3D "'V 3V3L>(1 - D 2 + D 4 ) 
2{D 2 -l) 2 (l + D 2 + D A ) , (1 + D 2 )u , 



3D(1 + L> 6 ) 



(! + £>> „ 

COS! 1) 



-{XP + PX) t=0 



(28) 



where 




To simplify things, we assume that the noise is a small effect compared to the free dynamics, i.e. 

r 



Using this approximation we get 



(Ho)(t) — huiexp 



Tlu 2 



((P 2 ) {t =o) + (X 2 ) {t =o)) 



(29) 



(30) 



(31) 



Again this result is equivalent to the one given by Savard et al. fici[| in the limit of white noise. In Fig. [j] both the exact 
(solid line) and the approximated (dashed line) solution are plotted for the vibrational frequency lj — 11.2(27r) kHz and 
the value of Tlo/2 is 0.1 in this plot. We choose the initial values to be (X 2 ) t=0 = {P 2 )t=v = (l/2)(XP+PX) t =o = 1/4, 
so that (HQ)t=o/TiLd = 1/2. 

IV. EFFECTS OF FLUCTUATING EXTERNAL ELECTRICAL FIELDS DURING QUANTUM GATE 

OPERATIONS 



We have investigated the overall effects of the fluctuations in the spring constant and the position in ion traps. It is 
of interest to look at the effects of those fluctuations on gate operations used for quantum computation jlij . The gate 
operations are performed by shining an additional laser, with a specific frequency and for a well determined time, on 
a two-level transition in an ion. This laser causes an interaction between the internal electronic states of the ion and 
the CM motion of the ion or of all the ions if there is more than one ion in the trap. 

We will denote the electronic qubit as \g), |e) for the ground state and the excited state respectively. The coding 
is such that the ground state is logical while the excited state is logical 1. The vibrational state will be denoted by 
the energy eigenstates 10)^, which are the ground state and the first excited state respectively. 

A controlled-NOT gate can be broken down into the circuit shown in figure |^, where a controlled phase shift is 
sandwiched between two ir/2 pulses with different phases on the target qubit, which in this case is the electronic 
qubit. The it/2 pulses produce rotations of the electronic qubit: 



U 



R ■ 



u 



R ■ 



|0)-1/V2(|0>-|1» 
|1>-1/V2(|0> + |1» 

|0>j-1A/2(|0> + |1» 
|l) i _-l/^(|0)-|l)) 



(32) 
(33) 



The controlled phase shift Up acts to produce a it phase shift only if both the electronic qubit and the vibrational 
qubit is in the logical state The total transformation from input to output is then given by 



\Vout) =U R UpU+p m ) . 
We take the most general input state to the controlled-NOT gate as 



(34) 



(a|ff>+/?|e»®(*|0) t 



ell) 



(35) 
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with a, 8, 7 and 6 being complex amplitudes satisfying 

M 2 + |/?| 2 = H 2 + H 2 = i. 

The first it/2 pulse Ut acting on the target (electronic) qubit, then produces the state; 

a 



l*i> 



S\g)\0) v + ^-=-5|e) 0)„ 



V2 



V2 



a + 8 8 



(36) 



(37) 



where the subscript 1 indicates that this is the state after the first rotation. The controlled phase shift operation just 
changes the sign of the last term. Then the final rotation takes this state to the output state 



\*out) = aS\g)\0) v + 36\e)\0) v + ae\e)\l) v + 8e\g)\l) v 



(38) 



If the vibrational qubit is subject to noise in the trap center position, i.e. there are fluctuating electrical fields, the 
gate will not operate as required. The noise does not effect the rotations, which only involve the electronic qubit 
(except for heating up the vibarional state during those rotations), however, it will disrupt the controlled phase shift 
operation which couples the electronic and vibrational systems. We first make the usual transformation into the 
interaction picture with 



Uq = exp 



-H Q t 



where 

Hq = Twja) a + TiUAcr+a- , 

to give 

Hi=H G - h\£(t) (a)e iu}t + ae - *"*) , 
where H G is the interaction picture Hamiltonian causing the gate operation and where 

A = (2hmujy 1/2 . 



(39) 



(40) 



(41) 



(42) 



To include the noise we need to calculate the effective controlled phase shift operation including the noise term over 
the time of the gate operation. The delta-correlated nature of white noise enables a simple approach in which noise 
terms can be separated and treated perturbatively while the gate interaction is treated to all orders. To enable this 
approximation we first transform to an interaction picture defined by 

!*(*)) =exp(-~H G t\ \*i(t)) . (43) 

The total time evolution for the controlled phase-shift gate over the gate time T is given by — U' P \^>i) with 

fT 



exp 



H G T T exp 



dt'H„ 



,(t') 



= U P U N [£(t)] 



where 



(t) = exp ( l -H G t (-ftA£(t) (aV" 4 + ae~ lut )) exp - l -H G t 



(44) 



(45) 



and where T is the time ordering operator and T the time required for the gate operation in the absence of noise. We 
have indicated the functional dependence on a particular realization of the noise by [£(£)]. The time-ordered evolution 
operator appearing as the second factor in Eq. ( |44| ) may be treated perturbatively in the stochastic amplitude, £(t), 
by carrying the Dyson expansion to second order. 

There are two ways to view the effects of the noise. One way is at the level of a single realization of a gate operation. 
This view enables us to see what the error states will be in the presence of noise. The second way is to determine the 
result of a gate operation by ensemble averaging the noise. This view enables us to give an average fidelity for the 
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gate operation in the presence of the noise. The two views will be referred to as a quantum trajectory picture and 
an ensemble average picture, respectively. In the quantum trajectory picture the output state is a functional of the 
particular noise history, £(t), over the gate operation time T, The output state is then a pure state of the form 

mm]) = UnU P U N [t(t)]\*i) > ( 46 ) 

where, as earlier, is the state after the first rotation of the electronic qubit. To find the likely error states the 
noise operator can now be expanded in powers of the noise amplitude. 

In the ensemble picture we need to average over all possible realizations of the noise over the gate time T . The 
output state is now a mixed state given by 

Pout = U R U P p[U P U R (47) 

and 

p'i = [ u N m)\*i)(*i\uim)p[tm > (48) 



where -P[£]d[£] is the probability functional for each noise realization. To calculate p[ we expand the unitary operator 
i/jv [£(*)] to second order in the noise and then average over the classical stochastic variables. The evolution of the 
density operator over the gate time T is then given by a Dyson-expansion |13| which to second order is 



p[=pi + ^J [H nolse {ti),p 1 }dt 1 + (J^j J dt x dt 2 [H noise (t 1 ), [H noise (t 2 ),pi]] ■ (49) 

Since E(dW(t)) = 0, the average over the second term vanishes and we only have to calculate the average over the 
third term. This is done by noting that fig] 



rt 



G(t')dW(t') / G(t")dW(t") = / G{t'fdt' . (50) 



We can quantify the effect of noise on the average through the fidelity defined by 

f(7) = (^out\p' out \^out) , (51) 

where \^ ut) is the output state for a noiseless gate operation while p' out is the output state of the gate averaged over 
all realizations of the noise. The fidelity is the probability that the system is in the desired state and will depend on 
the noise correlation strength 7. Substituting Eqs. (]34|, E?j) we find that 

F(j) = (tfiKltfi) • (52) 

To proceed we take two examples for realizing a controlled-NOT gate between the vibrational state and the internal 
state for one ion in the trap. 



A. Mutual phase-shift gate 

In the discussion of the perfect controlled-NOT gate we saw that the essential two-qubit operation is a controlled 
phase shift. We will discuss two different ways by which this can be done. The first way involves a mutual conditional 
phase shift of the vibrational and electronic degrees of freedom. This operator commutes with the vibrational quantum 
number. The second way involves an auxiliary electronic level and is used in the NIST scheme to produce a controlled- 
NOT gate fli"|| . We first discuss the mutual phase shift gate as this is simpler. 

The mutual phase shift gate is defined by the unitary transformation 

Up = cxp (— wra'a <£> |e) (e|) . (53) 

To include the noise we need to calculate the effective controlled phase shift operation including the noise term over 
the time of the gate operation, where the gate Hamiltonian is 

H G = hnata® |e)(e| (54) 
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and with the gate operation time such that kT = w, where k is a constant. An interaction of this kind can be produced 
by a carrier frequency excitation of the ion fl5| . 

The total time evolution is then given by Eq. ( flij ) with 

Hnoise — — 7iA£(i) (a exp [— int\e}(e\ — iu>t] + h.c.) , (55) 

where h.c. indicates the hermitian conjugate of the prccccding term. The effect of the noise is determined by expanding 
the second time ordered factor in Eq. (M4f) to second order, thus 

U N = l + i{av{T) + <Jv\T))-\ f dt\a < ^ + a^ ( ^P-\{au{t)+a^\t)} , (56) 



2 J | dt dt 
where the operators v are defined in the following way: 



v(t) = A^7 / tW(t')e~* (a,+re|e><e|) *' . (57) 



We first determine the likely error states in the quantum trajectory picture. The desired gate operation, Eq. (p33|), 
commutes with the phonon number operator a^a and cannot change the vibrational quantum number, and thus the 
vibrational states always remains in the logical basis of |0), |1). The noise factor Un, however, is linear in a and a 1 
and thus does change the phonon number. Keeping terms to second order in the stochastic amplitude means that, 
in a single realization of the gate, the effect of the noise is to 'leak' coherence into the vibrational states |2)„, |3)„ as 
well as changing the weighting of the vibrational qubit basis states \g), |e). For example, in the quantum trajectory 
picture, a single realization of the noise takes the state Eq. ( |37|) to the state 



V2 



K| 5 )|0), + V2^| 3 >|2),) + ^e^|e)|0) u + v / 2^| e )|2) t ,) + ...) , (58) 



where we have suppressed the second order terms, and the subscript 2 indicates that this is the state after the 
controlled phase shift gate operation, and the prime indicates a state corrupted by noise. The first term is the correct 
output state and all subsequent terms correspond to an error. The error terms are multiplied by random variables 
defined by 



T 



v g = \Vl dW(t)e- wt (59) 
Jo 

Ve = AV7 / dW{t)e- %{uJ+K ^ . (60) 



If the gate time T is lar ge c ompared to the vibrational frequency, these random variables have zero mean and the 
correlation functions are |16|| 

E(u g u g ) = E(u e u e ) = (61) 
E{v* gVg ) = E{v* e v e ) = A 2 7 T (62) 

E(y v e ) = . (63) 

As mentioned above we can calculate the fidelity of the gate operation, Eq. ( |5l| p2[ ) . The result up to first in 7 is 
F(T K ) = 1 - 2r K (l + 2|e| 2 ) + r„(l - |e| 2 )|e| 2 



|/3 - a| 4 1 + sin — cos(— + 2A) 

\ 7T(IjJ + K) K K J 

I \ K . t^Tt , . LU7T . \ .„ .,, „, An . L07T . . . L07T „ A . 

\a + 0\*[l+ — sin — cos — + 2A - P~a 2 a + /3 2 — -cos — sin — + 2A 

' ' V 7TUJ K K ) S 11 ' 7r(2w + K) K k' V K ' 



(64) 



where 



K hmuK 
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is the now dimensionless noise parameter and A is the phase difference between S and e 

A = (f> 5 - £ , (66) 

where we denote the phase of 8 (e) by (f>s (</> e ), respectively. 

As expected for an expansion up to second order in ^Fj the fidelity depends linearly on 7, since we have averaged 
out the first order term. The inverse dependance on the nonlinear coupling constant k is easy to understand. If this 
parameter is large, the gate time can be made very short in which case the noise has less time to act and produce 
an error. Thus the fidelity should approach one as k becomes large. Note that in general the fidelity depends on the 
initial state. The dependence of F on the initial state is plotted in Figure [|. The plot parameters are u> — 11.0(27r) 
MHz, Sir) = 1.0(27r) kHz, A = 0, and T a = 0.02. 



B. NIST Gate 



The Hamiltonian required to describe the controlled phase shift operation in the NIST gate is 

^G-?i^(|e)(aux| a t + |aux>(e|a) (67) 

where 77 is the Lamb-Dicke-parameter, and a and a? are the creation and annihilation operators for the quantized CM 
motion of the ion. This Hamiltonian is used to describe a 27r blue-sideband pulse between the excited state of the 
electronic qubit (logical 1) and an auxiliary state |aux). The rotations of the electronic states are carried out by 7r/2 
pulses that are assumed to act without noise. This assumption is reasonable since the effective Rabi frequencies for 
those transitions are higher than those for the blue sideband pulses, thus leading to much shorter pulse durations. 

As before we separate the pure gate operation from the noise by transforming to an interaction picture through the 
gate Hamiltonian. We then determine the average fidelity for the gate operation by averaging over the noise. 

Calculating p[ using Eq. (|49| ) and in particular Eq. ( f45|) is a rather tedious, but straightforward process. The fidelity 
rate turns out to be 



F(T a ) =l-2r„ [l + |e| 2 ] -r o |e| 



where 



-r a | e | 2 (l-| £ | 2 )l/3-a| 4 sin(^ 



AlTUJ 

~~Shi~ 



■2A 



Sir) 



Sir, 



Sir) 



+r a | e | 2 (l-| £ | 2 )|a + ^ 



fir? / 47rw \ 

1 H - sin — — 

27TCJ \ Sir) J 



_8ttw 167r(f2?7 - 2uo) 16n(Slr) + 2uj) 



Sir) 



2A 



+T a \e\ 2 (l-\e\ 2 )\a + f3\ 2 \f3- 



2 / 4iru \ ( 47rw 

sin I I COS I 



a sin 



\ Sir) J \ Sir) 

47T7 



2A 



Sit) 



Sir) 



n(Slr) - 4jjj) n(Slr] + 4u>) 



T a = 



hmtvSlr) 



(68) 



(69) 



is now the new dimensionless noise parameter and A is again the phase difference between the two vibrational states 
as defined above in Eq. (|6"6|). 

Again the fidelity depends linearly on r a . The dependence of F on the initial state we want to perform the gate 
on is plotted in Fig. ^. The chosen paramters for the plot are u> = 11.0(27r) MHz, Sir) = 12.0(27r) kHz, A = 0, and 
r„ = 0.02. 



V. DISCUSSION 



We have determined the heating rates due to a fluctuating trap potential and due to fluctuating electrical fields for 
the mena motional energy of an ion in a rf Paul trap. The potential use of ion traps as simple quantum computers 
in view we have calculated the effects of fluctuating electrical fields (considered one of the major sources of noise at 
the moment) during a controlled-NOT gate operation. We derived fidelities for two different ways of performing the 
required conditional phase shift needed for those gates. This analysis is particularly useful for the application of ion 
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traps to quantum computation. It gives an estimate on how strong those fluctuating fields can be to still perform a 
computation with a certain accuracy. 

Taking the current heating rate of the COM mode for the NIST ion trap 0, which is about 19 phonons per ms 
and assuming that these heating rate is due to fluctuating electrical fields (the reasons for those heating rates are 
not clear yet, so we just assume at this stage until experimentalists will come up with more elaborate data), we get 
a rough estimate for r a ss 0.02. With that value we get fidelities above 90 % for one gate operation, which certainly 
needs improvement to allow for more than one gate operation. 
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FIGURES 



FIG. 1. Plot of the exact (solid line) and approximated (dashed line) solution for the increase of the phonon number in 
time. This plot is for a vibrational frequency to — 11.2(27r) kHz and for Flo/2 — 0.1. We choose the initial values to be 
(X 2 ) t = = (P 2 )t=„ = (1/2) (XP + PX) t=0 = 1/4, so that <tf )t=o/(M = 1/2- 

FIG. 2. Schmeatic representation of a controlled-NOT gate: The controlled phase shift Up is sandwiched between the two 
single qubit rotations U£ and . 



FIG. 3. Plot of the dependence of F on the initial state we want to perform the gate on. The chosen parameters for the plot 
are w = 11.0(2tt) MHz, k = 1.0 MHz, A = 0, and F K = 0.02. 



FIG. 4. Plot of the dependence of F on the initial state we want to perform the gate on for the NIST gate. The chosen 
parameters for the plot are u> = 11.0(2tt) MHz, Qr] = 12.0(2tt) kHz, A = 0, and T a = 0.02. 
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